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Abstract. - Superslow diffusion, i.e., the long-time diffusion of particles whose mean-square 
displacement (variance) grows slower than any power of time, is studied in the framework of the 
decoupled continuous-time random walk model. We show that this behavior of the variance oc- 
curs when the complementary cumulative distribution function of waiting times is asymptotically 
described by a slowly varying function. In this case, we derive a general representation of the laws 
of superslow diffusion for both biased and unbiased versions of the model and, to illustrate the 
obtained results, consider two particular classes of waiting-time distributions. 
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Introduction. — Random processes in physical, bio- 
logical, social, economic and other systems often exhibit 
the anomalous diffusion behavior in the sense that at long 
times the variance cr 2 (t) of these processes increases non- 
linearly with time [TH3]. Usually, the variance follows the 
power law, i.e., o~ 2 {t) oc i Q , and two types of anomalous 
diffusion, subdiffusion (when < a < 1) and superdiffu- 
sion (when a > 1), are distinguished. The class of systems 
with such behavior of the variance is vast and growing. 
Subdiffusion has been observed, e.g., in amorphous solids 
[4], percolation clusters [5] and living cells [HE], and su- 
perdiffusion in turbulent flows 019], optical lattices [10] 
and foraging animals [TTHT2"] . 

In general, the time-dependence of the variance is not 
restricted by a power function. For example, under certain 
conditions the power law can be modified by a logarithmic 
function of time (see, e.g., refs. [T3HT5] ). But since the log- 
arithmic function varies much more slowly than the power 
one, this modification is not so important. In contrast, in 
the case of superslow diffusion, which is defined by the 
condition that a 2 (t)/t a — > as t — > oo for all a > 0, the 
variance does not follow the power law at long times. The 
Sinai diffusion for which a 2 (t) oc ln 4 i [16] represents the 
best known example of this type of diffusion. Some other 
examples of superslow diffusion have been found in resistor 
networks [17) . continuous-time random walks (CTRWs) 
[18], charged polymers [19], aperiodic environments [20] . 
iterated maps [5T], Langevin dynamics [52], fractional ki- 
netics [53], etc. A common feature of all these examples 



is that the long-time dependence of the variance, i.e., the 
law of superslow diffusion, is given by a power function of 
the logarithm of time: cr 2 {t) oc ln I/ t (v > 0). 

A natural question is if there exist laws of superslow 
diffusion that differ from a 2 {t) oc In" t. Since the defining 
condition a 2 (t)/t a — > holds not only for this particular 
law, one can expect that other laws also exist. In this 
letter we use the decoupled CTRW model to show that if 
the waiting-time distributions vary slowly at infinity (see 
below), then the corresponding laws of superslow diffusion 
form a wide class of slowly varying functions that increase 
without limit as time evolves. 

Description of the model. - Within the CTRW 
formalism (see, e.g., refs. [IH2]), the position X{t) (A(0) = 
0) of a diffusing particle is described by the continuous- 
time jump process 
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where N(t) S (0, 1, 2, . . .) is the random number of jumps 
up to time t (if N(t) = then X(t) = 0) and x n G 
(— oo, oo) is the random magnitude of the n-th jump. The 
number of jumps N(i) is characterized by the waiting 
times T n 6 [0, oo), i.e., random times between the particle 
jumps. These times and the jump magnitudes Xji are as- 
sumed to be independent and identically distributed with 
the probability densities p(r) and w(x), respectively. In 
the decoupled CTRW model it is also assumed that two 
sets of variables r n and independent of each other, 
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so that the joint probability density of these variables is 
given by the product p{r)w(x). 

In Fourier-Laplace space, the probability density P(x, t) 
of the process X(t) satisfies, in the decoupled case, the 
Montroll- Weiss equation [24] 
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where the Fourier and Laplace transforms are defined as 
Uk = J-{u(x)} = dxe lkx u(x) and v s = £{v(t)} = 
J °° dte~ st v(t) (Res > 0), respectively. Using eq. ©, the 
Laplace transform of the m-th (m = 1,2,...) moment of 
X{t) can be written in the form 



(X m (t)) s = (-»)' 
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with the angular brackets denoting an average over the 
sample paths of the random process X(t). Applying to 
eq. ([3]) the inverse Laplace transform, v(t) = £~ 1 {v s } = 
(l/27rc) f C _*°° dse st v s (it is assumed that the real param- 
eter c exceeds the real parts of all singularities of v s ), one 
obtains (X m (t)) = C- 1 {{X m (t)) s }. 

In this work we focus on the long-time behavior of the 
variance of the particle position, defined as 



\t) = (x*(t)) - (x(t)y 



(4) 



under the additional condition that all fractional mo- 
ments of the waiting-time distribution are infinite, i.e., 
T p = J °° drT p p(r) — oo at p > 0. By analogy with [25], we 
call the probability densities p(r) satisfying this condition 
the super-heavy-tailed densities. It should be noted that 
a special class of these densities with p(r) ~ aj(r ln 1+I/ r) 
(a > 0, t — > oo) has already been put forward by Havlin 
and Weiss [T5]. But here we consider a much wider class 
of super-heavy-tailed probability densities that are char- 
acterized by the asymptotic behavior 



p(r) ( T 
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where a positive function h(r) is slowly varying at infinity 
in the sense that the condition h(pr) ~ h(r) holds for all 
p > as r — > oo. Since the waiting-time distributions 
are assumed to be normalized, 7o = Jq°° drp{r) — 1, the 
admissible functions h(r) form a subclass of slowly varying 
functions that at r —> oo tend to zero fast enough. 

Laws of superslow diffusion. — In order to ensure 
the finiteness of the variance er 2 (i), we assume that the 
first two moments of w(x), l\ — J _ dxxw(x) and I2 — 
J^ oo dxx 2 w(x), exist. In this case, using the definition of 
the Fourier transform of w(x), one obtains 



w k = l+ihk- -l 2 k 2 + o{k 2 ) (fc ->•()), (6) 



and so eq. ([3]) for m = 1 and m = 2 gives 
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(X 2 (t)) s = 



hPs 



s(l-p s ) 2 s(l-p s ) 
(7) 

Our next step is to find the asymptotic expressions for 
(X(t)) s and (X 2 (t)) s when a positive parameter s ap- 
proaches zero. The motivation is that if these expressions 
satisfy certain conditions then the long-time behavior of 
the moments (X(t)} and {X 2 (t)) can be found directly 
from the Tauberian theorem for the Laplace transform. 
Let us first represent the Laplace transform of p{r) as 
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Ps = l- dr(l - e- ST )p(r). 
Jo 

Then, introducing the probability 

V(t) = / d T p(r) 



(8) 



(9) 



that the waiting time exceeds t and using the new variable 
of integration q = sr, we can reduce eq. ([8]) to the form 



Ps = 1 - / dqe~ q V 
Jo 



(10) 



The function V(t), known as the complementary cumula- 
tive distribution function of waiting times (or exceedance 
probability), satisfies the initial condition 1^(0) = 1 and 
tends to zero in a slowly varying way as t — > 00. The 
last property of V(t), which is crucial for our approach, 
can be straightforwardly verified using the asymptotic for- 
mula ([5]) and the condition h(pr) ~ h(r) (r — > 00): 
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(t — > 00). We note in this connection that the waiting- 
time cumulative distribution function F(t) = f Q drp(r) = 
1 — V(t) also varies slowly at infinity. 

Since V(pt) ~ V(t) as t — > 00, the function V(q/s) in 
eq. (fl"U|) at s — > can be replaced by V(l/s), yielding 



p,~i-y(-) (s^o). 



(12) 



Utilizing this result in eq. ([7]) , we find, in the limit of small 
s, that 
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if li = 0. Next, to find the asymptotic behavior of the 
moments (X(t)) and (X 2 (t)), we apply the Tauberian the- 
orem for the Laplace transform, which states (see, e.g., 
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ref. [26 ) that if the function v(t) is ultimately monotone 
and v s ~ s-TL(l/s) as s -> then u(t) - r^L^/r^) 
as t — > oo. Here, L(t) is a slowly varying function at in- 
finity, r(7) is the gamma function, and < 7 < 00. In 
accordance with this theorem, eqs. (fTB")) and (|14l) lead to 
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and 



(*(*)) = 0, (X 2 (t)) 
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(15) 



(16) 



respectively. We emphasize that in the biased case, when 
l\ 7^ 0, the main terms of the asymptotic expansion of 
(X(t)) 2 and (X 2 (t)) are different, see eq. (JT5]). The same 
situation occurs also for the biased subdiffusion, while 
for the biased normal diffusion (X 2 (t)) ~ (X(t)) 2 and 
(X 2 (t)) - (X(t)) 2 oc t [T31|27]. These features of the mo- 
ments arise from the fact that, since in the biased case 
P-ks 7^ Pks, the left and right tails of the probability den- 
sity P(x, t) arc different. 

Thus, using eqs. (gj), (H"5l) and (|T6"|). for the long-time 
behavior of the variance one immediately obtains 



a 2 (t) 



with 



V(t) 



ll/V 2 (t), h^Q 
k/V(t), h = 



dT-h(r). 



(17) 
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The asymptotic formula (|17l) , which describes both biased 
(h ^ 0) and unbiased (l\ = 0) diffusion in systems with 
super-heavy-tailed distributions of waiting times, is the 
central result of this letter. The most important is the 
fact that the exceedance probability V(t) varies slowly 
at infinity. It is this feature of V(t) that is responsible 
for the superslow character of diffusion in these systems. 
Indeed, it is well known [28] that if some function, say 
V(t), varies slowly, so does any power of this function. 
Hence, the function <r 2 (t) is also slowly varying, and thus 
the condition a 2 (t)/t a — > holds for all a > as t — > 00 
[28 . In other words, the variance cr 2 {t) grows to infinity 
(cr 2 (t) — s> 00 because V(t) —> as t — > 00) slower than 
any positive power of t, i.e., superslow diffusion occurs. It 
is worth noting that the slower the exceedance probabil- 
ity decreases, the slower this diffusion. At the same time, 
since a 2 {t)\i 1= £ / a 2 {t)\i 1= o — > 00 as t — » 00, the biased 
superslow diffusion is much faster than the unbiased one. 

Because the exceedance probability V(t) is slowly vary- 
ing, the Karamata representation theorem [28j suggests 
that V(t) at t — > 00 can always be written in the form 



V(t) ~ wexp( — j dr— e(r) 



(19) 



Here, v is a positive parameter and n is a non-negative pa- 
rameter which, without loss of generality, can be taken to 



be zero. In general, for an arbitrary slowly varying func- 
tion V(t) it is only required that the function e(r) tends to 
zero as t — y 00. But in our case V(t) — > in the long-time 
limit. Therefore, the function e(r) must be positive and 
satisfy the conditions e(r) — > and J dr'e(r')/r' — > 00 
as t — y 00. The asymptotic formula (|17[) with so defined 
exceedance probability V(t) represents the most general 
form of the diffusion laws in the reference model. If the 
asymptotic behavior of the function h(r) is known, then 
the asymptotic behavior of the function e(r) and the mag- 
nitude of the parameter v can be determined as 



e(r) 



h(r) 
V(r) 



(t 00) 



and 



lim exp 



dr'-e(r'))V(r) 



(20) 
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with V(t) taken from eq. (|T8|) . It should be emphasized 
that because the functions h{r) and V(r) are slowly vary- 
ing, their ratio, the function e(r), also varies slowly. This 
property holds for all waiting-time probability densities 
characterized by the asymptotic behavior §5§ . 

We complete this section by noting that the laws of 
superslow diffusion (|17[) are obtained within the simplest 
CTRW model. The fact that these laws include the dif- 
fusion law a 2 (t) oc In" t as a very special case does not 
mean, of course, that this model describes in detail the 
systems mentioned in the Introduction. In this context 
we recall that our goal is not the application of the de- 
coupled CTRW model to those systems but rather the use 
of this model to show that the laws of superslow diffusion 
constitute a broad class of slowly varying functions. 

Illustrative examples. — As a first example of super- 
heavy-tailed distributions of waiting times whose asymp- 
totic behavior is described by eq. ([5]), we consider a two- 
paramctric class of such distributions characterized by the 
probability density 



(r — 1) ln r 1 rj 
(?y + r) ln r (77 + r) 



(22) 



with the parameters r > 1 and 77 > 1 (these conditions 
guarantee that p(r) is positive and normalized). In this 
case the exceedance probability ([9]) is calculated exactly, 



V(t) = 



( ln?7 



VMi7 + *) 



(23) 



and thus V(t) ~ (In 77/ lni)*" -1 as t — > 00. The substitu- 
tion of the last asymptotic result into eq. (jTT)) leads to the 
following asymptotic formula: 



a 2 (t) 



(/i/ln r - 1 7 ? ) 2 ln 2(r_1) i, h^O 



h =0, 



(24) 



which for l\ — was first obtained in ref. 18.. Finally, 
since according to eq. (|2"2"|) h(r) ~ (r — 1) ln r_1 77/ In" r as 
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t — )• oo, from eqs. (|2T))) and (|2ip we find e(r) ~ (r— 1)/ In t 
and w = ln r ~ 77. 

Our second example is characterized by the waiting- 
time probability density of the form 



P (t) = 



(r — l)(lnln?7) r 



(77 + t) 111(77 + t) [In ln(r? + r)] r 



(r > 1, r\ > e). The main feature of this density is that its 
right tail is heavier than in the previous case, i.e., p{r) at 
r — > 00 tends to zero slower than the probability density 
from eq. (l2"2"j) . Therefore, the exceedance probability 



V(t) 



In In 77 
lnln(rj + t) 



(26) 



which corresponds to the probability density (|25|) . de- 
creases slower than V{t) from eq. f[23]) and, as a conse- 
quence, the variance (|17j) increases with time slower than 
in the previous case: 

_ f [^/(lnlnr/r-^^lnlnt) 2 ^- 1 ), 1^0 

a 2 it) ~ <^ 27 
\ [; 2 /(lnlnr;)''- 1 ](lnlnt)'- 1 , h=0. 

Taking also into account that according to eq. (|25|) h(r) ~ 
(r — l)(lnln77) r_1 /[lnr(lnlnr) r ] as r — > 00, for this case 
we obtain e(r) ~ (r — l)/(lnrlnlnr) and tj = (lnln?7) , '~ 1 . 

The above examples are illustrative and, of course, the 
diffusion laws (j2"4")l and (|2T|) do not exhaust all possible 
laws of superslow diffusion. As is shown in this work, 
the possible laws form a class of slowly varying at infin- 
ity functions cr 2 (t) satisfying the condition a 2 (t) -> 00 as 
t — > 00. This class is sufficiently large because it contains 
infinitely many functions which at t — > 00 grow both faster 
and slower than a given function from this class. For ex- 
ample, the function a 2 (t) oc exp(h/ t)j hit (0 < j3 < 1) 
grows faster than tr 2 (t) from eq. (|24|) . and the function 
a 2 (t) oc In In In t grows slower than a 2 (t) from eq. (|27|) . It 
should also be stressed that, although for the probability 
densities p(r) characterized by the condition ([5]) the mean 
waiting time is infinite, the regime of superslow diffusion 
is reached at finite times exceeding the longest time scale 
of the corresponding density. In particular, the diffusion 
laws ([Ml) an d (EH) take place at t » 77. 



Conclusions. — Using the decoupled continuous-time 
random walk model characterized by the waiting-time dis- 
tribution functions that vary slowly at infinity, we have de- 
termined a wide class of the superslow diffusion laws. This 
class consists of the slowly varying functions, describing 
the variance of the particle position, which tend to infin- 
ity slower than any power of time. We have established the 
connection between the long-time behavior of the variance 
and the asymptotic behavior of the exceedance probabil- 
ity. Specifically, in the cases of biased and unbiased su- 
perslow diffusion the variance is inversely proportional to 
the second and first powers of the exceedance probability, 
respectively. Due to this difference, the biased superslow 



diffusion at long times is much faster than the unbiased 
one. Finally, by applying the Karamata representation 
theorem, we have found the most general form of the laws 
of superslow diffusion in the reference model. 
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